Tunable Superfluidity and Quantum Magnetism with Ultracold Polar Molecules 
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By selecting two dressed rotational states of ultracold polar molecules in an optical lattice, we 
obtain a highly tunable generalization of the t-J model, which we refer to as the t-J-V-W model. In 
addition to XXZ spin exchange, the model features density-density interactions and novel density- 
spin interactions; all interactions are dipolar. We show that full control of all interaction parameters 
in both magnitude and sign can be achieved independently of each other and of the tunneling. As 
a first step towards demonstrating the potential of the system, we apply the density matrix renor- 
malization group method (DMRG) to obtain the ID phase diagram of the simplest experimentally 
realizable case. Specifically, we show that the tunability and the long-range nature of the interac- 
tions in the t-J-V-W model enable enhanced superfluidity. Finally, we show that Bloch oscillations 
in a tilted lattice can be used to probe the phase diagram experimentally. 

PACS numbers: 67.85.-d, 71.10.Fd, 33.80.-b, 71.10.Pm 



Experiments with ultracold atoms have recently ex- 
tended the range of candidate systems for realizing un- 
conventional states of matter and enabled the simulation 
of models describing condensed matter phenomena [1]. 
One major goal of current research at this interface be- 
tween condensed matter and atomic physics is to emulate 
the Heisenberg and t-J models, which are believed to un- 
derlie certain quantum magnetic materials [2] and high- 
temperature superconductors [3], respectively. However, 
in the ultracold atom realization of these models, the 
small superexchange interaction J [1] makes the under- 
lying physics extremely challenging to observe. At the 
same time, ultracold polar molecules [4], such as KRb 
[5, 6] and LiCs [7], have recently been produced in their 
electronic and rovibrational ground states. In this Let- 
ter, we show that when such molecules are localized in 
an optical lattice, their rotational degree of freedom can 
be used to simulate tunable Heisenberg-like models at 
unit filling of the lattice and, in the presence of dop- 
ing, a fully tunable generalization of the t-J model that 
we refer to as the anisotropic t-J-V-W model. Dipole- 
dipole interactions that give rise to this model are orders 
of magnitude stronger than superexchange interactions 
J in ultracold atoms and can therefore better compete 
with other relevant energy and time scales such as, for 
example, those responsible for decoherence. Moreover, 
we show that the resulting long-range interactions are 
fully controllable with DC electric and continuous-wave 
microwave fields. As a first step towards demonstrating 
the potential of the model, we use DMRG [8] to obtain 
the ID phase diagram for the simplest experimentally rel- 
evant case and show that, at low fillings, the superfluid 
phase is enhanced relative to the one in the conventional 
t-J model [9]. We propose to probe the phase diagram 
using center-of-mass Bloch oscillations. Given that KRb 



has already been loaded into a 3D lattice, our proposal 
is applicable to current experiments. 

The Hamiltonian and its features. — We consider di- 
atomic polar molecules in their electronic and vibrational 
ground state partially polarized by a DC electric field 
along z, confined to the x-y plane [10], and loaded in 
that plane into the lowest band of a square optical lattice. 
As described below, microwave fields are used to isolate 
in each molecule two dressed rotational states |too) and 
|mi) and to obtain the t-J-V-W Hamiltonian 
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The two terms describe tunneling and dipole-dipole inter- 
actions, respectively; () denotes nearest-neighbor bonds. 
Specifically, ct creates a fermionic [5] or bosonic [6, 7] 
molecule on site j (position R 3 ) in dressed rotor state 
m € {mo,77ii}; in this Letter, we focus on fermions. 
Large reaction rates [10] between two molecules on the 
same site enforce the hardcore constraint. We define 

n jm = c j m C jrrn n j = ^2 m n jmi Sj' = c jm C i'"'l' = 

( n jm a ~ n jm 1 )/'^ and use units in which h = 1. The J z , 
V, and W terms can be understood by thinking of \mo) 
and \mi) as classical permanent dipoles oriented along 
z, while the \mo)-\mi) transition dipole moment gives 
rise to the J± term. We tune these dipole moments by 
constructing dressed states |m ) and \m\) out of bare ro- 
tor states using microwave fields. This, in turn, allows 
for the full controllability of J z , J±, V, and W, which 
is one of the main results of the present Letter. Lattice 
Hamiltonians based on more than one molecular rota- 
tional state have been considered before in Refs. [11-22]. 
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FIG. 1. (color online), (a) Eigenenergies of Ho — BN 2 — doE 
as a function of E (top axis uses d and B of KRb). Only 
AI > states are labeled. (b,c) Examples of level con- 
figurations {\m )_, |mi)}: (b) _{Vo|0) + VI - a|T>, (c) 
{Va|2) + \/T^|2), Vb|l) + Vc|T)+Vl - 6 - c|2)>. Red (blue) 
levels make up the dressed rotor state |mo) (|mi)). 



An important difference of Eq. (1) from the Hamiltonian 
studied in Ref. [17], which is most closely related to our 
work, is the presence of the J± term [11]. 

Eq. (1) possesses several aspects that can lead to un- 
conventional many-body phases. Most notably, all the 
interactions are long-range; in particular, repulsive long- 
range density-density interactions are expected to stabi- 
lize superfluid correlations [23]. Furthermore, the novel 
W term can break SU(2) symmetry of the system even 
when J z — J±. It is also crucial that J z and J± can be 
tuned in both sign and magnitude (up to ~ 100 kHz in 
LiCs) independently from each other and from t. This 
contrasts with the cold atom realization of the t-J model, 
where J <C t. Finally, we note that, in the limit of unit 
filling (n = 1), only the terms with J z and Jj_ survive, 
realizing an XXZ-model with long-range dipolar interac- 
tions, and, for J z — J±, the Heisenberg model. 

Tuning J z , Jj_, V, and W . — We now show how to use 
DC electric and microwave fields [12, 14-18, 22, 24-28] 
to control J z , J±, V, and W, leaving a detailed deriva- 
tion of Eq. (1) to Ref. [29]. A polar molecule in a DC 
electric field is described by the rigid-rotor Hamiltonian 
Hq = £>N 2 — doE, with the rotational constant B, the 
angular momentum operator N, and the dipole moment 
operator d. We set d p — e p ■ d, where e = z and 
e± = t(x ± iy)/V% [25]. At E = 0, H has eigen- 
states | A, M) obeying N 2 | A, M) = N(N+1)\N, M) and 
N Z \N,M) = M\N,M). As E is increased, eigenstates 
with the same M mix, forming superpositions of \N, M) 
with different N . We denote the state at finite E adia- 
batically connected to \N,M) as \4>n,m) and define, for 
notational convenience, \N) = \4>j? t o), \N) = |0jv,i)j an d 
\N) = \<f>N,2) [see Fig. 1(a)]. For appropriate choices 
of levels, energy conservation ensures that dipole-dipole 
interactions between two molecules preserve their total 
N z , resulting in the many-body interaction Hamiltonian 

#dd = | Eiyy iH^-RjI-^^^ + ICdWd^+dW^)]. 
We now focus on the two configurations of levels and mi- 
crowave fields shown in Figs. l(b,c), which allow, respec- 



tively, for the realization of the simple case J z — V — 
W = and for the full tunability of J z , J±, V, and W. 

In the first configuration, states \mo) — \/a\0) + 
\/l — a\l) and \mi) = |1) [Fig. 1(b)] are chosen, where 
the coefficient a in the dressed state |mo) is controlled 
by the ratio between the Rabi frequency and the de- 
tuning of a (7 + -polarized microwave field acting on the 
|0)-|1) transition [26]. In 40 K 87 Rb, microwave inten- 
sity of a few W/cm 2 is sufficient to address all hyper- 
fine levels equally. Projecting on states |0), |1), and 
|1) and keeping only energy conserving terms, we ob- 
tain d^d ( Q j) ={il 1 \0l)(l0\ i j+h.c. + ® k J2 s Vs\s)(s\k and 

d W d W) +d w d w = -^^ioi^xoj.. _ M 2_| 1I)(I1 | y . + h c 5 

where k £ {hj}> s G {0,1,1}, Ms = (s|do| s ), an d 
Mss' = ( s \d p \s'} for the appropriate p. The minus sign 
in front of /u, 2 ^ and /u 2 ^- is crucial to the tunability of V, 
W , J Zl and J± and appears because two dipoles rotat- 
ing in the x-y plane interact on average attractively [26] . 
Projecting on |mo) and |mi) and again keeping only en- 
ergy conserving terms, we obtain 
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where p,q G {0, 1} and A = a/i + (1 — a)p,j, A± = fii, 
B a = -MqT ^ _ a )' Bl = 0, Jj_ = 2/i^a - /i 2 T (l - a). 
A p can be thought of as the effective dipole moment of 
\tn p ), while B p is the contribution to Hdd from tran- 
sition dipole moments between rotor states making up 
\m p ). From the comparison of this equation and Eq. 
(1), we can read off V = \{A Q + Ai) 2 + B + B x ]/A, 
W = {Al + B Q -A\-B x )/2, J z = (Ao-A^ + Bq + Bl 
A simple case J z — W = V — and J± > studied 
below can be approximately implemented using a = 1, 
which does not require a microwave field. A small E 
field (dE/B > 0.1 in 40 K 87 Rb) is needed to prevent 
dipole-dipole and hyperfine interactions from populating 
|0i, ±i). At dE/B = 0.1, V, W, and J z are two orders of 
magnitude smaller than J± and can be neglected. 

The second configuration we consider is {|mo), |^i)} = 
{>|2) + VT^|2), Vb\l) + v^|l) + y/l-b-c\2)} [Fig. 
1(c)]. The three microwave fields shown in the fig- 
ure, which allow to control the coefficients a, b, and 
c, together with the DC electric field constitute four 
knobs that allow for the full control over the four co- 
efficients V, W, J z , and J± [30]. We find that in a small 
sphere in the 4-dimensional (dE/B,a,b,c) space around 
the point (dE/B, a, b, c) = (2.97, 0.059, 0.56, 0.38), where 
J z = J± = V = W = 0, one can achieve any value of 
V, W, J z , and J± up to an overall positive prefactor. 
Similarly, in the special case where a single microwave 
couples |1) and |1), we have a two-dimensional (dE/B, b) 
subspace (with a — 0, c = 1 — b), in which J z and J± can 
be fully controlled for simulations of the XXZ model at 
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FIG. 2. (Color online). Phase diagrams of the t-J± chain with 

(a) nearest-neighbor and (b) dipolar Jj_. We identify a metal- 
lic (repulsive Luttinger liquid) phase with dominant spin- 
density- wave correlations (SDW), a gapless superfiuid with 
dominant triplet and singlet superfiuid correlations (TS/SS), 
a singlet superfiuid with a spin gap (SS+SG), and phase sep- 
aration (PS). The spin gap is < 0.35i in (a) and < 0.7t in 

(b) . Solid lines indicate phase transitions (green: SG closes; 
red: inverse compressibility becomes zero); dashed lines are 
extrapolations. The numbers show the value of K p on the 
dotted lines. The line K p = 2 is a crossover line within 
the SS+SG phase [31]. The shaded region in (b) displays 
K p = 1 ± 0.15 as an estimate of the numerical accuracy. 



unit filling - around (dE/B, b) — (4.36, 0.56), where they 
both vanish. While these examples prove full controlla- 
bility in their respective cases, for any desired relation- 
ship between V, W, J z , and J±, there is likely a different 
level configuration that gives stronger interactions and 
uses weaker E, lower microwave intensity, and/or more 
convenient microwave frequencies. 

Phase diagrams of the nearest-neighbor and dipolar t- 
J± chains. — The full tunability of the t-J-V-W model 
provides access to a great variety of models with poten- 
tially exotic physics. As the simplest example of this 
physics, we present in Fig. 2 the ID phase diagram in 
the limit V — W — J z = 0, which is one of the simplest 
experimentally achievable cases (see above). Before an- 
alyzing dipolar interactions, we present in Fig. 2(a) the 
phase diagram of the nearest-neighbor t-J± chain as ob- 
tained using DMRG and following the analysis of Ref. [9] . 
The diagram is qualitatively similar to that of the stan- 
dard t-J chain [9]: At fillings n < 0.65, we identify a 
repulsive Luttinger liquid (Luttinger parameter K p < 1) 
with dominant spin-density-wave (SDW) correlations, an 
attractive Luttinger liquid [K p > 1) with dominant sin- 
glet and triplet superfiuid correlations, a singlet super- 
fluid with a spin gap, and phase separation. At larger 
fillings, the spin gap is always zero, but the other phases 



In Fig. 2(b), we present the phase diagram of this 
system in the presence of dipolar interactions. At low 
fillings, SDW, the gapped singlet superfiuid, and phase 
separation are obtained. This suggests that experiments 
with ultracold molecules can be a useful tool for explor- 
ing the phase diagram of the standard t-J model for ar- 
bitrary values of J in contrast to J <C t in ultracold atom 
realizations. At the same time, crucially, both diagrams 
in Fig. 2 feature a significant enhancement of the super- 
fluid region compared to the original t-J model because 
the absence of attractive density-density interactions sup- 
presses phase separation. Furthermore, the maximum 
value of the spin gap in the dipolar t-J± chain is twice 
that in the nearest-neighbor t-J± chain, which is, in turn, 
twice that in the original t-J model. The larger spin gap 
should facilitate its experimental observation. Note that, 
as a conservative estimate of the numerical errors in com- 
puting K p in the presence of long-range interactions with 
up to 100 sites, we estimate the true line K p — 1 to lie in 
the shaded region 0.85 < K p < 1.15. Thus, the gapless 
superfiuid cannot be identified in Fig. 2(b) within our 
numerical precision. Furthermore, while the line K p = 1 
and the line where the spin gap closes may coincide, our 
analysis of the correlation functions cannot rule out the 
existence of an exotic intermediate phase with a spin gap, 
K p < 1, and dominant superfiuid correlations. 

Preparation and detection. — Ground-states at specific 
points in the phase diagram can be prepared, e.g., by 
applying an additional microwave field coupling |m ) and 
|mi) and performing an adiabatic passage from an easily 
accessible state to the desired ground state by tuning the 
Rabi frequency and the detuning of the microwave field 
[18]. Direct probing of molecules [32] and conversion of 
molecules back to atoms [5] can in principle both be used 
for detection via noise-correlations in the time-of-flight [1] 
or via in-situ single-site imaging [33, 34]. 

As a specific example of a detection technique avail- 
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FIG. 3. (color online). Bloch oscillations in a dipolar t-J± 
chain of 20 sites with open boundary conditions at fillings 
n — 0.1 and 0.2, and tilting field E t \\ t /t — 1 per site. The main 
plot shows the difference between the center-of-mass position 
xcm at time zero and at the first minimum. The vertical line 
indicates the approximate value of J±/t, at which the spin 
gap closes for n = 0.2. The inset shows the time evolution of 
iom at n = 0.2 at three indicated values of J±/t. 



4 



able in current experiments, we propose to use Bloch 
oscillations [35, 36]. In Fig. 3, we present our results 
as obtained via the Krylov-space variant of the adaptive 
t-DMRG [37, 38] for 20 sites when adding a linear field 
along the chain at filling n = 0.1 and n = 0.2 for different 
values of Jj_/t. In the singlet supcrfluid, the amplitude 
of the oscillations drops and the frequency of the oscil- 
lations doubles (see J±/t = 7 in the inset) relative to 
the gapless phase due to the presence of bound pairs. 
Bloch oscillations should be observable in direct absorp- 
tion imaging [32]. While neither the frequency nor the 
amplitude of the oscillations show any sharp features at 
this small system size, the fit of experimental data to nu- 
merical results should allow for the location of the phase 
transition (vertical line in Fig. 3) even for small system 
sizes. As a complementary method for identifying the 
transition, we propose spectroscopic measurement of the 
spin gap [39]. 

Outlook. — We have presented a toolbox for simulating 
a highly tunable anisotropic t-J-V-W model with polar 
molecules. The advantages of this molecular toolbox over 
its atomic counterpart are higher energy scales and in- 
dependent tunability of interactions and tunneling. This 
toolbox should enable the simulation of condensed matter 
phenomena, as well as the stabilization and controllable 
preparation of unconventional phases, such as d-wave su- 
pcrfluids [3]. The phase diagram of the experimentally 
simple case of a t- J± chain shows an enhanced superfluid 
region, which we propose to probe via Bloch oscillations. 
We expect that the 2D t-J-V-W model can similarly be 
tuned into exhibiting enhanced superfluidity [40]. 

The present Letter also opens other exciting research 
avenues. In particular, natural extensions of the model 
include [29]: spatially anisotropic interactions produced 
by a tilt in the DC electric field, spin-dependent tunnel- 
ing obtained by adjusting lattice beams, S > 1/2 models 
realized by choosing more than two dressed states, and 
systems with an orbital degree of freedom encoded in the 
nuclear spin. Moreover, by considering molecular Wigner 
crystals [24, 41], where the intcrmolccular distances are 
smaller than in an optical lattice, one can further increase 
the interaction strength. Furthermore, by analogy with 
Ref. [18], we expect our ideas to be extendable to Ry- 
dberg atoms. Finally, one can envision applications of 
the present system to quantum computation (especially 
if one uses nuclear spin to store information), precision 
measurements, and controlled quantum chemistry [4]. 
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